function y = func_alfa_1(p, CONSTS)
    eps = CONSTS.eps;
    g = CONSTS.g;
    q1 = func_q1_from_p(p, CONSTS);
    
    y = -1 + (p.^2+q1.^2-eps)/g;
end
